The Nature of Quantum Hall States near the Charge Neutral Dirac Point in Graphene 
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We investigate the quantum Hall (QH) states near the charge neutral Dirac point of a high mobility 
graphene sample in high magnetic fields. We find that the QH states at filling factors v — ±1 depend 
only on the perpendicular component of the field with respect to the graphene plane, indicating them 
to be not spin-related. A non-linear magnetic field dependence of the activation energy gap at filling 
factor u = 1 suggests a many-body origin. We therefore propose that the u = and ±1 states arise 
from the lifting of the spin and sub-lattice degeneracy of the n — LL, respectively. 

PACS numbers: 73.63.-b, 73.21.-b, 73.43.-f 
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The experimental observation of the Quantum Hall 
(QH) effect in single atomic sheet of graphene [l|, 0] has 
attracted much attention recently, particularly due to 
the unique electronic transport observed in this mate- 
rial. Compared with the conventional integer QH effect 
in many other two-dimensional (2D) systems, the Hall 
resistance (Rxy) quantization condition in graphene is 
shifted by a half integer: R^y = ±gs{n + l/2)^, where n 
is the Landau level (LL) index, e is the electron charge, h 
is Planck's constant, and = 4 is the LL degeneracy, ac- 
counting for spin and sub-lattice symmetry in graphene. 
This quantization condition leads to the QH effect ap- 
pearing at filling factors i^ = ±2,±6,±10,---. It is now 
understood that this unique QH effect is related to the 
Quasi-relativistic nature of the charge carriers in graphene 
0i 0i J stemming from the unusual linear dispersion re- 
lation of its bands near the charge neutral Dirac point in 
the graphene band structure @,13|- 

More recently, the QH effect in graphene has been 
studied in the extreme quantum limit in a very strong 
magnetic field Q. New QH states, corresponding to fill- 
ing factors v = 0, ±1, ±4, are clearly resolved in magnetic 
fields S > 20 T, indicating a lifting of the four-fold degen- 
eracy of the n = LL and a two-fold degeneracy of the 
71 = ±1 LLs respectively. While angular dependent ac- 
tivation energy gap measurements indicate that the QH 
states at ly ~ ±4 are spin states, the origin of the QH 
states at ly = 0,±1 remains unresolved. 

The nature of these QH states near the charge neu- 
tral Dirac point is of fundamental interest. There have 
been numerous theoretical investigations [ O l. llCt I 111 . Il2l . 

H, [3, [H, M, [13, El, M, M, M, of these 

states, and their origin is currently under considerable 
debate. Recently, Abanin et al. [26[ suggested that the 
= QH state is spin-polarized and dissipative, owing 
to counter-propagating edge states at the charge neutral 
point, supported by a finite metallic resistivity at low 
temperatures. 



In this Letter, we present an experimental investiga- 
tion of the QH states near the Dirac point. We find that 
the ly ~ ±1 states depend only on the out of plane com- 
ponent of the applied magnetic field, and show a rather 
large energy gap with an approximately square root de- 
pendence on the magnetic field. This suggests a many- 
particle origin of this splitting as it would originate from 
the breaking of the sub-lattice degeneracy of the n = 
LL at the Dirac point. As a consequence, and by elimi- 
nation, our results would imply that the v = QH gap 
is induced by the lifting of the spin degeneracy. 

Our sample is a high quality graphene specimen with 
mobility as high as ~ 2 x lO''^ cm^/Vs measured at car- 
rier density rig = 4 x 10"'^^ cm~^. The graphene sheet 
is mechanically extracted from Kish graphite following 
a method similar to the one described in [H, The 
sample is deposited onto a Si substrate, which serves as 
a gate electrode separated from the sample by 300 nm 
of insulating Si02. To perform transport measurement, 
multiple electrodes are patterned in van der Pauw geom- 
etry (inset to Fig. 1) using conventional electron beam 
lithography, followed by Au/Cr (30/3 nm) thermal evap- 
oration and a standard lift-off process. The electronic 
transport is measured over the temperature range of 4.2- 
300 K, using a lock-in technique. The sample is mounted 
on a single-axis tilting stage to allow in situ tuning of 
the angle, 9 ~ cos^^ {Bp/ Btot), where Btot is the total 
magnetic field and Bp is the component perpendicular to 
the graphene plane. 

The four-fold degeneracy of the n — LL of graphene 
consists of a two- fold degeneracy from the spin symmetry 
and a two-fold degeneracy from the sub-lattice symmetry. 
One may be able to distinguish the origin of any partic- 
ular splitting by performing magnetotransport measure- 
ment in a tilted field, where a spin splitting depends on 
Btot, whereas a sub-lattice splitting (caused by electron- 
electron correlations) would only depend on Bp. Figure 
1 shows the measured magnetoresistance, Rxx, with re- 
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FIG. 1: (color online). Magnetoresistance, Rxx, with respect 
to the back gate voltage, Vg, at a temperature oi T = 4.2 
K and in two different magnetic fields: Btot ~ 45 T, Bp = 
20 T (dashed curve); and Btot = 30 T, Bp = 20 T (solid 
curve). The minimum magnetoresistance J?™" substantially 
increases as Btot decreases for u = —4 state, while for the 
v — ill states -RJi?" practically does not depend on Btot for 
the same 



spect to the back gate voltage, Vg, at a temperature of 
T = 4.2 K and in two different magnetic fields: Btot = 45 
T, Bp = 20 T (dashed curve); and Btot = 30 T, Bp ^ 20 
T (solid curve). Note that in this comparison we use the 
same Bp but different Btot- The minimum magnetoresis- 
tance i?™" substantially increases as Btot decreases for 
1/ = —4 state, while for the = ±1 states i?""" hardly 
varies with Btot for the same Bp. This observation sup- 
ports the spin related origin of the v = ±4 splittings [8| , 
and also, more importantly, suggests that the = ±1 
states are likely due to the breaking of the orbital degen- 
eracy of the sub-lattice symmetry in the n = LL. 

In order to further characterize the nature of ;/ = ±1 
QH states, we measure the activation energy of i?™" 
at fixed magnetic fields. Figure 2(b) displays the Ar- 
rhenius plots of R 



of the V — \ state [27|. A well- 
defined thermal activation behavior is readily observable 
(i?™™ ~ exp[-Ai;/2fcBT], where ks is the Bohzmann 
constant), and the corresponding energy gap, Ai?, can 
be extracted for different magnetic fields. In Fig. 2(a), 
we plot the obtained energy gaps dX v = 1, denoted as 
/S.E{v = 1), as a function of B-field. For comparison, 
we have also reproduced the measured AE{v = ±4) of 
the same sample from Ref. 0. We find /S.E{v = 1) to 
be considerably larger than lS.E{v = ±4). For instance 
at 45 T, A£'(i/ = 1) > Al:\E{v = ±4). Moreover, unlike 
for l^E{v = ±4), which showed a linear _B-field dcpen- 
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FIG. 2: (color online), (a) The measured activation energy 
gap, /S.E, as a function of magnetic field for the QH states 
at filling factors u = 1 (square), v = —4 (solid dot), and 
V = A (open dot). While AE{i' = ±4) are linear with respect 
to B-field (data are reproduced from Ref. [^), AE{v — 1) 
can be fitted better by a y/B dependence (solid curve), (b) 
Arrhenius plots of R^x minimum of the v = 1 state for three 
different magnetic fields: 45 T, 37 T, and 30 T. The straight 
lines are linear fits to the data. 



dence, the i3-field dependence of /S.E{u — 1) does not 
seem to follow such a simple relationship. Forcing a lin- 
ear fit onto AE{i' ~ 1) produces a positive y-intercept 
which would indicate an unphysical, negative LL energy 
width. In fact, a \/B behavior provides a better fit to 
the AE{i' = 1) data as shown in the solid curve in Fig. 
2(a). 

The lack of a linear dependence and the existence of a 
roughly \/B dependence point to a non-spin origin and 
possibly to a many-particle origin of the gap, and sug- 
gest that the = ±1 states are associated with a spon- 
taneous breaking of the sub-lattic e sy mmetry driven by 
the electron-electron interactions O 0, 0, 3, 13, E, 
11,0, EH. In this picture^£;(z/ = ±1) is expected to 
be on the scale of e^/ds [16[, and thus proportional to 
Vb, where Ib = \/ti/eB is the magnetic length and e is 
the dielectric constant. We calculate that e^/els ~ 1100 
K, assuming e = 4 and i? = 45 T. This value is much 
larger than the Zeeman energy Ez = gfJ-sB ~ 60 K at 
B ~ 45 T, where 5 = 2 is the ^-factor and /xb is the 
Bohr magneton. This simple evaluation indicates that 
/S.E{v = ±1) >> Ez in the experimentally accessible 
magnetic field, suggesting the importance of electron- 
electron interaction under mangetic fields. 

In Fig. 3, we summarize our current understanding 
of the sequence of the QH states near the charge neu- 
tral Dirac point of graphene in a schematic of the LL 
hierarchy. We use up-arrows and down-arrows to rep- 
resent the spin of the charge carries, and solid (blue) 
and open (red) dots for different valleys in the graphene 
band structure. Since our measurements suggest that the 
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FIG. 3: (color online). Schematic of the LL hierarchy in 
graphene in magnetic fields. The up and down arrows repre- 
sent the spin of the charge carries, and the solid (blue) and 
open (red) dots indicate different valleys in the graphene elec- 
tronic band. 



ly = ±1 states are associated with the valley splitting of 
the n = LL due to electron-electron correlations, the 
QH state at ly ~ must be related to the spin split- 
ting of this LL. However, we also note that the behavior 
of Rxx and R^y at = is completely different from 
that of any other QH states away from the charge neu- 
tral Dirac point. Unlike usual QH states, the i/ = QH 
state does not show a resistance minimum in R^x nor a 
clear resistance plateau in Rxy This state only becomes 
visible as a plateau in the Hall conductance. Figure 4 
displays Rxx vs. Vg near the Dirac point over a wide 
range of temperatures. No activation behavior has been 
observed at = 0. Recently, Abanin et al. provided a 
possible interpretation for the existence of this state as 
the consequence of counter-circulating edge states with 
opposite spin (26j . Such a state would be consistent with 
our proposal shown in Fig. 3. 

We now address the relative size of the energy gap 
between the levels displayed in Fig. 3. In a sufficiently 
large magnetic field, typically > 20 T, the QH states 
in high-quality graphene specimens are robust, even at 
room temperature [28[ (see also Fig. 4). From the inset 
to Fig. 4, we estimate the activation energy gaps of the 
QH states at v = ±2 and B = 45 T to be AEiv = 2) « 
890 K (solid dots), and AE{iy =^ -2) w 570 K (open 
dots). In a simplistic view, we may interpret these values 
as the gap between the n = and the n = ±1 LLs in 
Fig. 3. In a single particle picture with un-lifted spin and 
sub-lattice degeneracies [29|, the LL energy spectrum in 
graphene can be described by 



E„ = sgn(7i) 



2ehv\B \n\ 



(1) 



where vp \s the Fermi velocity of graphene with a typical 
value of VF ~ 10*^ m/s [1 S, 0, Ell. Hence, at B = 45 
T, the calculated energy spacing between the n = and 
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FIG. 4: (color online). Magnetoresistance, Rxx, as a func- 
tion of the back gate voltage, Ki, at B = 45 T in a wide 
temperature range, from room temperature to liquid helium 
temperature. A robust QH effect is clearly seen even at room 
temperature. No activated behavior has been observed in Rxx 
of the = state (at Dirac point). Inset: Arrhenius plots of 
Rxx minimum of the v = ±2 states at B = 45 T. The straight 
lines are linear fits to the data. 



the n = 1 LL would be Eq^i « 2800 K. This value is 
more than three times larger than the largest measured 
energy gap AE{i' = 2) w 890 K. We believe that such 
a large discrepancy cannot be accounted for in terms of 
spin or sub-lattice symmetry splitting, nor by any rea- 
sonable LL broadening. We first rule out the possibility 
of an enhanced spin splitting since the g-factor is not en- 
hanced by exchange in a completely filled LL. With the 
bare value of g = 2, the spin splitting of the LL reaches 
only gusB ^ 60 K at 45 T. We further eliminate the sub- 
lattice degeneracy splitting as a potential explanation of 
the observed discrepancy, as this gap AE{v = ±1) must 
collapse when the Fermi energy lies between the n = 
and the n = 1 LL. Finally, the LL broadening due to 
scattering may lead to a reduction of the energy gap. 
However, we infer such a reduction is negligible compared 
to Eq^i ~ 2800 K: from our low temperature measure- 
ments (T = 1.4 K) we estimate a LL broadening of 
r « 20 K. Since the mobility of graphene changes only 
by ^ 30% from 30 mK to room temperature 32l|. F is 
likely irrelevant on the scale of Eq^i- 

At this stage we do not know what causes this strong 
reduction of the measured LL energy gap as compared 
to the calculated one. One may speculate that many- 
particle effects are partially responsible for the discrep- 
ancy. Recent infrared experiments in graphene 
35[ between LL levels ti = and n = ±1 yield a rather 
good agreement with Eq. 1. Yet theory suggests that 
as much as 30% of this energy is due to many-particle 
corrections [s^. Since both single-particle energy level 
and many-particle correction have -v/B-ficld dependences 
in graphene, they cannot be separated. This would sug- 
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gest that the Fermi velocity (typically vp ~ 10^ m/s) 
can be as much as ~ 30% smaller in reality, the dif- 
ference being made up for by many-particle corrections. 
If activation energy measurements were much less af- 
fected by such corrections, the expected gaps would fol- 
low the bare value of vp, leading to a ^ 30% reduction 
from Eq. 1. This would bring the calculated energy 
gap of i5o-ti ~ 1900 K closer to the measured gap of 
/S.E{v = 2) « 890 K, but still leaves a substantial, unre- 
solved discrepancy. The implication of such a reinterpre- 
tation oi vp would be considerable and more extensive 
studies of the LL spectrum in graphene will be required 
to verify such a trend. 

In conclusion, we study the QH states in graphene 
at filling factors = ±1 in tilted magnetic fields and 
elevated temperatures. Our results indicate that the 
V = ±.1 QH states originate from the lifting of the sub- 
lattice symmetry of the n = LL caused by electron- 
electron interactions. Measurements of the activation en- 
ergy gaps of the QH states near the Dirac point indicate 
a significant deviation from a simplistic single-particle 
model, which suggests that many-particle effects need to 
be taken into account to understand the LLs near the 
charge neutral Dirac point. 
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